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Three-dimensional topological insulators (TIs) are characterized by spin-polarized 
Dirac-cone surface states that are protected from backscattering by time-reversal sym¬ 
metry. Control of the spin polarization of topological surface states (TSSs) using 
femtosecond light pulses opens novel perspectives for the generation and manipula¬ 
tion of dissipationless surface spin currents on ultrafast timescales. Using time-, spin-, 
and angle-resolved spectroscopy, we directly monitor for the first time the ultrafast 
response of the spin polarization of photoexcited TSSs to circularly-polarized fem¬ 
tosecond pulses of infrared light. We achieve all-optical switching of the transient 
out-of-plane spin polarization, which relaxes in about 1.2 ps. Our observations estab¬ 
lish the feasibility of ultrafast optical control of spin-polarized Dirac fermions in TIs 
and pave the way for novel optospintronic applications at ultimate speeds. 


The emerging field of ultrafast spintronics in 
condensed-matter physics relies on the possibility of 
achieving efficient control of pure spin-currents, spin- 
polarized electrical currents and spin-configurations on 
ultrafast timescales m- One alternative for this pur¬ 
pose is the use of technologically relevant information- 
storage devices which are composed by ferromagnetic 
layers, and utilize laser-assisted switching for ultrafast 
remagnetization. For instance, the use of circularly- 
polarized femtosecond (fs) laser pulses has been estab¬ 
lished as a promising route to excite and coherently con¬ 
trol spin dynamics in magnets without involving spin pre¬ 
cession or external magnetic fields [1]. The efficiency 
of such magnetic devices might be even enhanced by 
manipulation of isolated spins with pulses as short as 
the timescale of the exchange interaction, and within a 
single-photon shot [1]. 

On the other hand, the growing demand of low-power 
consumption by using the least possible current in such 
devices, has motivated in recent years a completely dif¬ 
ferent avenue to achieve generation and control of pure 
spin currents at ultimate speeds. This alternate pathway 
relies on the possibility of controlling carrier spins in low¬ 
dimensional systems through the spin-orbit interaction 
[MSI- The process of spin-current generation and its ma¬ 
nipulation on sub-picosecond (ps) timescales is based on 
non-magnetic materials, and control of the electron spin 
is achieved solely by optical means mu- This method 
also allows for the generation of spin-polarized electrical 
currents where in addition to the spin there is net flow of 
charge. With the advent of new classes of electronic ma¬ 
terials such as topological insulators (TIs) [TUHT^ . this 
unique route towards ultrafast optical control of spin cur¬ 
rents and spin-polarized electrical currents appears to 
be very promising as it might lead to much lower en¬ 
ergy consumption as compared to devices entirely com¬ 
posed by ferromagnetic layers. The key difference re¬ 
sides on the fact that, while insulating in the bulk due to 
strong spin-orbit coupling, the surface electronic struc¬ 


ture of TIs is characterized by spin-polarized Dirac-cone 
topological surface states (TSSs) that are protected by 
time-reversal symmetry, hence by unprecedented prop¬ 
erties such as forbidden backscattering [TfUITK] . This 
phenomenon opens the way for the realization of dis¬ 
sipationless spin-current devices or ultrafast information 
transport based on spin-polarized electrical currents on 
the surface of these materials [7]. For instance, it has 
been recently shown that such spin-polarized electrical 
currents flowing across TI surfaces can be utilized to ex¬ 
ert spin-transfer torque effects in adjacent ferromagnetic 
layers, being the size of the torque greater than the one 
induced by any other material so far |16j . 

In this context, generation of ultrafast spin- 
polarization transport in TIs means the creation of a 
nonequilibrium spin population of charge carriers above 
the Fermi level in sub-ps timescales [IIlIIH]. The unique 
response of photoexcited TSSs to fs light pulses has been 
demonstrated in a number of different experiments [ 13 - 
[2H], primarily using time- and angle-resolved photoe¬ 
mission spectroscopy (tr-ARPES) [TUH?7] . Based on a 
pump-probe scheme providing access to transient pop¬ 
ulation and depopulation dynamics, tr-ARPES has be¬ 
come the most powerful tool in systematically revealing 
the nonequilibrium dynamics of TI materials in energy 
and momentum space. Recent studies utilizing this tech¬ 
nique have primarily identified bulk-assisted electron- 
phonon scattering as the underlying mechanism respon¬ 
sible for the ultrafast momentum relaxation of photoex¬ 
cited carriers in TIs m [lo]. It has been shown that 
following optical excitation by intense infrared fs pulses, 
this mechanism leads to an ultrafast response of both 
bulk and surface electronic structure that is character¬ 
ized by coherent-phonon oscillations at terahertz frequen¬ 
cies |25) . In addition, exotic topological quantum phases 
such as transient Floquet-Bloch states have been experi¬ 
mentally observed, and coherent control of these photon- 
dressed surface bands has been achieved by circularly- 
polarized ultrashort midinfrared pulses with energies be- 
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low the bulk band gap [22]. 

Although extensive spin-resolved ARPES studies on 
various TI compounds have successfully identified the he¬ 
lical spin texture of TSSs naniiMsi], or the condi¬ 
tions for its manipulation using circularly-polarized pho¬ 
tons in equilibrium [33RS5| . the critical response of the 
spin properties of photoexcited TSSs to the incident light 
and its polarizations on ultrafast timescales has remained 
unexplored, mainly due to the lack of spin resolution in 
previous time-resolved experiments. Based on indirect 
methods such as time-resolved reflectivity or Kerr rota¬ 
tion, it has been suggested that relaxation of spin and 
charge degrees of freedom of photoexcited TSSs occurs on 
very different timescales m- The experiments indicate 
that spin relaxes much faster than charge, in particular 
within the first few fs following optical excitation, which 
corresponds to the timescale of electron-electron scatter¬ 
ing [MMl [3SI- Such an effect enabled the observation 
of light-induced electrical currents flowing across TI sur¬ 
faces by means of transport measurements, and without 
the need of time resolution [18]. Specifically, reversal of 
these currents with the helicity of circularly-polarized fs 
pulses allowed to attribute them to spin-polarized elec¬ 
trical currents originating from the unique helical spin 
texture of TSSs [T5]. However, a direct microscopic ob¬ 
servation of the ultrafast spin-polarization dynamics of 
TSSs shedding light on the time evolution and underly¬ 
ing origin of these helicity-dependent photocurrents, is 
still missing. 

Here, for the first time, we provide the key experi¬ 
mental evidence that allows us to directly determine the 
ultrafast nonequilibrium dynamics of the momentum- 
space spin configurations induced on TI surfaces using 
circularly-polarized light. We perform spin-resolved 
tr-ARPES measurements on the naturally p-doped TI 
Sb 2 Te 3 , which is characterized by completely unoccupied 
TSS bands in equilibrium. We demonstrate all-optical 
control of the nonequilibrium spin-polarization of Dirac 
fermions on ultrafast timescales and further identify the 
underlying mechanism that triggers the spin-relaxation 
process following optical excitation. 

Nonequilibrium bulk and surface dynamics 

To study the temporal evolution of the transient popula¬ 
tions in the unoccupied band structure of Sb 2 Te 3 , we per¬ 
form pump-probe experiments under the two-color opti¬ 
cal excitation scheme shown in Fig. la. Infrared and 
ultraviolet fs-laser pulses of 1.5 and 6 eV photon energy 
play the role of pump and probe beams, respectively, 
both incident on the sample at an angle of (j) = 45° with 
respect to the surface normal. The time resolution of 
the present experiment is ^200 fs, and the pump fluence 
~300 /iJ/cm^. We optically excite electrons from the 
occupied states into unoccupied bands near the Fermi 
level using circularly-polarized pump pulses of positive or 
negative helicity (C-|- or C-). Subsequently, we monitor 
the ultrafast response of the transiently modified elec¬ 
tron distribution by emitting photoelectrons above the 
vacuum barrier using linearly s-polarized probe pulses at 
variable time delays At. 

Figure Ib shows low-energy tr-ARPES measurements 


recorded along the TK direction of the Sb 2 Te 3 surface 
Brillouin zone (SBZ). The spectra are taken with this 
direction oriented perpendicular to the light incidence 
yz plane in Fig. la, and with infrared-pump pulses 
of C-l- polarization. At negative delays (At=-lps), we 
probe the electronic structure in equilibrium, as evi¬ 
denced by the lack of photoemission intensity above 
the Fermi level. Following optical excitation, two tran¬ 
sient electron populations in the initially unoccupied TSS 
and bulk-conduction band (BCB) are clearly observed. 
These populations rapidly thermalize to a nonequilib¬ 
rium Fermi-Dirac distribution by means of electron- 
electron scattering processes [IMD- Subsequently, the 
electronic system undergoes an ultrafast relaxation and 
completely disappears through an avalanche of hot elec¬ 
trons within ~3 ps (right panel in Fig. lb). This process 
is characterized by a significant energy transfer to the 
phonon system, which occurs via electron-phonon scat¬ 
tering m HOI 1211131] ■ 

Such a rapid electron recombination evolves mainly via 
two cooperative mechanisms, one being the effective ac¬ 
tion of the BCB population, which behaves as an electron 
reservoir continuously feeding the TSS bands nn ng, 
the other one an electron bottleneck effect arising at the 
Dirac node, which is located at an energy of Ed ~0.15 
eV. Differently from previous findings m I11I2H, the 
bottleneck effect is surprisingly strong in our data, and 
manifests as a more pronounced time-dependent accumu¬ 
lation of intensity above the Dirac node, near the bottom 
of the upper Dirac cone in Fig. lb. In contrast, the spec¬ 
tra further reveal a significant suppression of the Dirac- 
node spectral weight that is independent of At (see Sup¬ 
plementary). These effects can also be distinguished in 
measurements performed over a larger momentum region 
around the SBZ center (Figs. Ic and Id). Here we show 
a full-ARPES mapping of the TSS measured at At=500 
fs near the Dirac node (Fig. Ic), together with the cor¬ 
responding constant-energy surfaces extracted in a wider 
energy range (Fig. Id). The bottleneck effect emerges 
from a substantially reduced phase space for electron re¬ 
laxation at the Dirac node. Such a blockade of the hot- 
electron decay is expected to be considerably strong in 
the presence of a band gap at the Dirac node, especially 
if the gap exceeds the maximum phonon energy avail¬ 
able for electron relaxation (~25 meV). At first glance, 
this scenario might seem consistent with our data if we 
consider that, unlike linearly-polarized light, the rotating 
electric field of the circularly-polarized pump pulse might 
create a nonequilibrium massive-fermion state through 
the transient breaking of time-reversal symmetry [22] . 
However, for this condition to be fulfilled, it is established 
that a strong coherent coupling between the intense laser 
field and the TSS is primarily required, implying that 
the pump-pulse photon energy must not exceed the size 
of the bulk band gap (~0.25-0.3 eV) [52]. This limita¬ 
tion explains the absence of a dynamical gap at the Dirac 
node in Fig. lb, a conclusion which is further supported 
by the absence of a bright Dirac node in recent tr-ARPES 
experiments using linearly-polarized 1.5 eV photons [19] . 
These findings motivate us to attribute the Dirac-node 
suppression in our data to final-state effects caused by 
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FIG. 1: Ultrafast dynamics of Dirac fermions following optical excitation by circularly-polarized infrared pulses. 

(a) Experimental geometry, (b) tr-ARPES spectra obtained along the FK direction of the SBZ using pump pulses of positive 
helicity (C+). (c, d) Full-ARPES mapping 500 fs after photoexcitation. The Dirac cone exhibits circular constant-energy 
contours that evolve into warped regions with increasing energy, (e) Observation of direct optical transitions from deeper-lying 
bulk states into Dirac bands above the Fermi level. The TSS bands right below and in the immediate vicinity of the Fermi 
level strongly overlap with bulk-valence band states, (f, g) Normalized tr-ARPES intensities of (f) the TSS and (g) BCB 
as a function of pump-probe delay. The curves are extracted from the small energy-momentum windows shown in (e). (h) 
Energy-dependence of the TSS and BCB decay times for opposite pump-pulse helicities. (i) Transient electronic temperatures. 


matrix-elements introduced by the probe pulse. 

In this context, rather than coherent interaction with 
the TSS, our data reveal signatures of strong coupling 
between the circularly-polarized infrared pump pulse and 
bulk valence-band states located at an energy of ^ — 1 eV 
below the Fermi level (see Fig. le). This effect manifests 
itself as a pronounced intensity cut-off which remains at 
a constant energy of ^0.5 eV above the Fermi level up to 
At=400 fs, before electrons from higher levels relax down 
below this energy. The energy difference (vertical red- 
solid arrow in Fig. le), exactly corresponds to the pump 
photon energy, indicating that both the TSS and BCB 
are largely populated through direct optical transitions 
from these deeper-lying levels (black-dashed line in Fig. 
le). Besides the Fermi level widening, other transitions 
at off-normal wave vectors contribute with a small weight 
only (vertical red-dashed arrow in Fig. le). Such direct 
interband excitation generates a high concentration of 
valence-band holes which rapidly relax within the subse¬ 
quent dynamics. This process, which leads to an excess 
of hot electrons above the Fermi level, might represent 
another important contribution to the recently discov¬ 
ered electron-hole asymmetry in the Dirac cone |26j . 

Further insight into the dynamics of excited states can 
be obtained by analyzing the integrated tr-ARPES inten¬ 
sity within small energy-momentum windows (labelled 
from 1 to 17 in Fig. le). This method allows us to disen¬ 


tangle the TSS and BCB dynamics from the tr-ARPES 
spectra, as shown in Figs. If and Ig, respectively. As ex¬ 
pected, the hot-electron avalanche within the TSS bands 
gives rise to a delayed electron filling and progressively 
slower electron decay when the energy decreases (Fig. 
If). A similar but less pronounced behavior can be ob¬ 
served for the BCB population (Fig. Ig). Once the BCB 
population reaches its maximum, both the TSS and BCB 
dynamics are very similar, resulting in decay times that 
follow each other upon reversal of the circular light polar¬ 
ization (Fig. Ih). This effect pinpoints bulk-to-surface 
coupling as one of the driving forces underlying the ul¬ 
trafast electron relaxation. The strong coupling results 
in a slower electron thermalization and an average hot- 
electron temperature of ^2.3x10^ K (Fig. li). The 
rise times within the higher-energy windows are longer 
than the time resolution, indicating that multiple elas¬ 
tic events of interband electron-electron scattering addi¬ 
tionally play a role, in agreement with previous studies 

HHHMl- 

Up to At~600 fs, the electron decay is dominated by 
intraband scattering, which leads to a momentum re¬ 
laxation of the TSS and BCB populations from higher 
to lower energy, towards the bottom of each individ¬ 
ual band. The subsequent dynamics mainly depends on 
the effective electron transfer from the BCB to the TSS 
through interband scattering. At At=950 fs, due to such 
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a BCB-emptying process, a large portion of the BCB 
electrons has already decayed into the TSS bands (see 
Fig. lb). Both interband and intraband scattering pro¬ 
cesses are strongly mediated by the emission of acoustic 
and optical phonons lillTlIMl- 

Remarkably, the transient TSS population in Fig. lb 
is characterized by an asymmetric intensity distribution 
at opposite -l-fc||,x and wave vectors. The asymme¬ 
try originates from a circular-dichroic effect induced by 
the pump pulse polarization m- The dynamic recom¬ 
bination of excited BCB electrons within the Dirac cone 
leads to a progressive decay of this asymmetry with in¬ 
creasing time (see Supplementary), an effect that is the 
manifestation of a time-dependent electrical current or, 
in other words, a transient net flow of charge across the 
sample surface. The overall decay of the generated cur¬ 
rent is slightly less than the relaxation time of the total 
population, and we will return to this point later. At 
the onset of the optical excitation, such an imbalance in 
the density of excited carriers is equivalent to a circular- 
dichroic effect in the Dirac bands that reaches an absolute 
value of ^15% upon reversal of the pump-pulse polariza¬ 
tion from C-l- to C- (see Supplementary for more details). 
This effect, not present in the BCB, indicates that rever¬ 
sal of the pump-pulse helicity reverses the direction of the 
electric current, and that this phenomenon occurs only 
on the surface. These findings clearly demonstrate pre¬ 
vious conclusions from transport experiments m and 
differ from tr-ARPES studies using 1.5 eV photons, so 
far conducted with linearly-polarized pump pulses [191- 
miiiiiisiiis]. In this respect, we point out that a fully- 
symmetric transient TSS population with respect to the 
surface normal, such as the one generated using linearly- 
polarized infrared light m. would not give rise to the 
ultrafast electrical current observed here. The reason is 
that an entire suppression of the net flow of charge on the 
surface occurs if, for all values of At, there is an equal 
number of excited electrons moving in opposite momen¬ 
tum directions. 

The overall relaxation dynamics observed here pro¬ 
ceeds in a much faster timescale than the one previously 
found in other prototypical TIs such as Bi 2 Se 3 [19] and 
Bi 2 Te 3 |26j . where the nonequilibrium TSS population 
exhibited long lifetimes of more than 10 and 50 ps, 
respectively. We attribute this significant difference to 
the fact that, in the case of Sb 2 Te 3 , we observe a smaller 
separation in energy-momentum space between the TSS 
and BCB populations. This situation facilitates to a 
large extent phonon-mediated transfer of BCB electrons 
into the TSS bands, leading to a stronger bulk-to-surface 
coupling and thus to a pronounced bottleneck effect 
through an accelerated electron dynamics. 

All-optical sub-ps control of the spin polarization 

To investigate the transient response of the Dirac-cone 
spin polarization to circularly-polarized infrared light, 
we perform spin-resolved tr-ARPES measurements at a 
fixed time delay of 500 fs upon reversal of the pump- 
pulse helicity from C-|- to C- (Eig. 2). The tr-ARPES 
dispersion corresponding to this time delay is shown in 
Eig. 2a, where vertical-dashed lines indicate the elec¬ 


tron wave vectors at which spin-resolved measurements 
are presented. The spin-integrated energy-distribution 
curves (EDCs) extracted at these wave vectors are em¬ 
phasized in red in Eig. 2b. Each EDC contains a pro¬ 
nounced TSS peak located at an energy of ^0.33 eV, 
where the Dirac-cone distortion due to hexagonal warp¬ 
ing is rather weak (see Figs. Ic-ld). This peak is ac¬ 
companied by a less intense but clearly distinguishable 
high-energy tail due to the BCB contribution. 

Our spin-resolved measurements (Figs. 2c and 2d) 
reveal that circularly-polarized infrared photons tran¬ 
siently flip the photoexcited TSS spin polarization per¬ 
pendicular to the surface, and that reversing the pump- 
pulse helicity reverses the direction of the resulting out- 
of-plane spin polarization. These effects are clearly ob¬ 
served for a fixed wave vector in the in-plane and 
out-of-plane spin-resolved EDCs shown in Figs. 2cl-2c3, 
and in the corresponding net spin polarizations (Figs. 
2dl-2d3). Specifically, the in-plane component of the 
TSS spin polarization is nearly zero (Fig. 2dl), and the 
out-of-plane one reaches a maximum absolute value of 
^12% upon reversal of the pump-pulse helicity from C- 
to C-l- (Figs. 2d2 and 2d3). Moreover, the out-of-plane 
spin polarization does not reverse when going from —fey 
to -|-fc||,a; wave vectors (Figs. 2c4 and 2d4). These find¬ 
ings establish that the ultrafast net flow of charge on 
the surface is accompanied by a fully controllable spin 
polarization. 

In contrast, the transient BCB population contributes 
negligibly to the measured spin polarizations. This re¬ 
sult is consistent with previous conclusions from time- 
resolved reflectivity experiments using infrared light m, 
and with the fact that due to bulk-inversion symmetry, 
the BCB population is expected to be essentially unpo¬ 
larized. The effect is seen in our data as a nearly-zero 
spin polarization of the high-energy BCB tail above ~0.4 
eV. Although the pump-pulse helicity might initially cre¬ 
ate spin-polarized BCB carriers, it has been shown that 
due to strong spin-orbit coupling, their spin polarization 
will rapidly decay to zero faster than the time resolu¬ 
tion of the present experiment m- On the other hand, 
the impact of the spin-orbit interaction on the TSS is 
expected to cause spin relaxation rather than decoher¬ 
ence of its transient spin polarization. In particular, the 
theoretically expected spin-relaxation time of the TSS in 
prototypical TIs can exceed several ps [381ES]. The fact 
that the BCB in Fig. 2 is already unpolarized after 500 
fs indicates that in the presence of strong spin-orbit cou¬ 
pling, electron-electron scattering acts as the dominant 
depolarizing mechanism for bulk states. This finding im¬ 
plies that the electron-electron scattering rates for bulk 
and surface states are different by at least one order of 
magnitude. We note that this difference is not neces¬ 
sarily connected to an elevated hot-electron temperature 
of the BCB population. Moreover, differently from the 
TSS, the BCB spin polarization is prone to be affected by 
electron-impurity scattering events following optical ex¬ 
citation. Other effects associated with spatial diffusion 
of BCB electrons are likely to play a minor role mi¬ 
lt is difficult to directly compare the observed mag¬ 
nitude of the transient TSS spin polarization in Fig. 2 
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FIG. 2: Ultrafast spin-polarization control of Dirac fermions. Spin-resolved tr-ARPES spectra are obtained along the 
TK direction upon reversal of the pump-pulse helicity from positive (C -f) to negative (C-), at a fixed time delay of 500 fs. 
(a) Corresponding tr-ARPES dispersion indicating the momentum cuts at which spin-resolved measurements are presented 
(vertical-dashed lines), (b) Energy-distribution curves (EDCs) extracted from (a). EDCs corresponding to spin-resolved 
measurements are emphasized in red. (c) Spin-resolved spectra of the in-plane (cl) and out-of-plane (c2-c4) spin-polarization 
components obtained at opposite momenta using different light polarizations (C-|- and C-). The out-of-plane projection is 
perpendicular to the surface normal, and the in-plane one tangential to the circular constant-energy contour of the TSS. 
(dl-d4) Corresponding net spin polarizations. The in-plane components associated to panels (c3, c4) and (d3, d4) provide 
similar results as in panels (cl) and (dl). (e) Simplified schematics of the ultrafast optical orientation process, (f-h) One-step 
photoemission calculations revealing dominant transitions from s-like initial states to p-like final states. In (f) and (g), p- 
and s-like initial-state orbitals are suppressed in the transition matrix elements, respectively, (h) The normalized intensity 
difference reveals distinct contributions from p- (red) and s-like orbitals (blue). 


to the one expected from ground-state model calcula¬ 
tions |40[ I41j , as the spin polarization in the nonequilib¬ 
rium state can be strongly affected by multiple-scattering 
events that occur during electron thermalization and 
within the subsequent hot-electron dynamics, resulting 
in a substantial spin depolarization HZ]. For instance, 
the TSS spin polarization can be reduced within the 
timescale of electron-electron scattering and below the 
time resolution of the present experiment depending on 
the rate of electron-electron collisions [22]; and with¬ 
out affecting the transient charge current on the surface. 
Elastic events of interband electron-electron scattering 
giving rise to an enhanced bulk-surface coupling upon 
electron thermalization may additionally act as spin- 
depolarization channels m- However, in Fig. 2, we find 
a non-zero value of the transient TSS spin polarization, 
indicating that the fastest depolarization events associ¬ 
ated with elastic electron-electron scattering are not suf¬ 
ficient to completely suppress the spin polarization. The 
same is true for ultrafast spin-depolarization processes 
associated with an hyperfine interaction between excited 
electrons and nuclear spins (Sj , indicating that the hyper¬ 
fine coupling is too weak to transiently create an internal 
magnetic field. 

Another depolarization pathway is intraband decay of 
unpolarized electrons from energies higher than ~0.6 eV, 
where the TSS and BCB connect to each other bridging 
the gap of the volume |43| . This depolatization pathway 
has been proved very effective in recent spin-resolved tr- 


ARPES experiments focusing on the surface resonances 
of Bi 2 Se 3 using linearly-polarized light |44]. However, 
this contribution seems to be small in our data, as most 
of the excited electrons are generated below the intensity 
cut-off observed at ~0.5 eV in Figs. Ib-le. Moreover, 
spin depolarization of the Dirac cone can also occur upon 
momentum relaxation at every stage of the hot-electron 
avalanche due to scattering with phonons |9] . 

In Fig. 2e, we summarize the transient behavior of 
the Dirac-cone spin polarization by showing a simplified 
scheme of the optical spin orientation process. The circu¬ 
lar pump pulse aligns the hot-electron spin polarization 
antiparallel to the surface normal for C-|- and along the 
opposite direction for C- light, leading to a transient spin 
texture that is similar to the one recently predicted |33j 
and observed using circularly-polarized 6 eV photons in 
equilibrium [341135j . The difference in the present case is 
that the spin-reorientation process occurs for nonequilib¬ 
rium conditions and is solely caused by the pump-pulse 
helicity in sub-ps timescales. We emphasize that after 
the first photoemission step induced by the pump pulse, 
the linear polarization of the 6 eV probe pulse is perpen¬ 
dicular to the out-of-plane TSS spin polarization in the 
intermediate state, thus leaving the photoelectron spin 
polarization unchanged upon the second excitation above 
the vacuum level (see Supplementary for details). The 
spin-reorientation process induced by the pump pulse can 
be viewed in analogy to the circular-photogalvanic effect 
previously observed in semiconductor quantum wells |45j . 
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FIG. 3: Probing the spin-polarization and charge dynamics of Dirac fermions. (a,b) Ultrafast temporal evolution 
of the TSS (a) out-of-plane and (b) in-plane components of the spin polarization following photoexcitation with circularly- 
polarized infrared pulses of positive helicity (C-E). Similar to Fig. 2, the out-of-plane spin projection is perpendicular to the 
surface normal, and the in-plane one tangential to the circular constant-energy contour of the TSS. (c) The relaxation dynamics 
of the net out-of-plane spin polarization (blue) and charge asymmetry (red) exhibits different timescales. The spin polarization 
relaxes with a time constant of about 0.4 ps, while the charge asymmetry is delayed by ~1 ps. 


The process is the result of optical spin pumping into the 
unoccupied Dirac cone with electrons that are excited 
with preferential spin orientations. The resulting TSS 
out-of-plane spin polarization ultimately originates from 
the spin dependence of the relativistic dipole matrix ele¬ 
ments in the presence of strong spin-orbit coupling, which 
in turn depend on the symmetry of the states involved 
in the optical transition [35]. 

Hence, to examine the wave function character of the 
states that contribute to the spin-reorientation process, 
we perform one-step model photemission calculations un¬ 
der equilibrium conditions using linearly s-polarized 6 
eV photons (Figs. 2f-2h). We analyze the symmetry 
of the wave functions by suppressing p- or s-like initial- 
state orbitals in the corresponding transition matrix el¬ 
ements and normalize the intensity (Figs. 2f and 2g, re¬ 
spectively). While the TSS contains contributions from 
all p-orbitals as expected |40l ET], the calculations re¬ 
veal that the transient TSS spin polarization originates 
predominantly through direct optical transitions from 
spin-degenerate bulk states of s-like character (Fig. 2h). 
These states are thus represented by an initial-state wave 
function of high symmetry (angular momentum i=0), 
which is a precondition for optical spin orientation in 
solid surfaces using circularly-polarized light |3Sj. The 
calculations also reveal bulk-like surface resonances lo¬ 
cated at the BOB edges (black-dashed lines in Fig. 2h). 
We emphasize that the wave function of the TSS con¬ 
sists of a complex superposition of different spin direc¬ 
tions appearing on different atoms 03111 ]. As a con¬ 
sequence, the out-of-plane TSS spin polarization shown 
in Figs. 2d3 and 2d4 reveals only weak momentum de¬ 
pendence at opposite wave vectors, and preserves time- 
reversal symmetry. This process can be understood in 
terms of orbital angular-momentum transfer from the 
circularly-polarized pump pulse, which in the presence 


of strong spin-orbit coupling triggers spin-flip transitions 
as a function of the light helicity. Moreover, the magni¬ 
tude and direction of the spin polarization are expected 
to be time-dependent |38l [39] , implying in this way that 
transient changes in the relative orbital occupations are 
equally possible. 

Ultrafast spin-polarization dynamics 

To understand the spin-polarization dynamics of the 
TSS, in Fig. 3 we monitor the temporal evolution of 
the out-of-plane (Fig. 3a) and in-plane (Fig. 3b) com¬ 
ponents of the spin polarization after optical excitation 
with circularly-polarized infrared pulses of positive he¬ 
licity (C-E). Similar to Fig. 2, the measurements are 
taken at a fixed +fc||,a wave vector and along the FAT 
direction of the SBZ. We observe a rapid and progres¬ 
sive suppression of the out-of-plane spin polarization, 
while the in-plane component remains nearly zero. After 
^ 1.2 ps, the out-of-plane spin polarization is strongly re¬ 
duced, implying that spin decoherence due to scattering 
with phonons rather than a gradual recovery of the ex¬ 
pected Dirac-cone helical spin texture likely dominates 
the momentum-relaxation process. Moreover, the fact 
that the spin polarization disappears faster than the TSS 
population reveals that the transient spin-polarized elec¬ 
trical current originating from Dirac fermions transforms 
into a pure electrical one. 

In Fig. 3c, we gain further insight into this transfor¬ 
mation by comparing the out-of-plane spin-polarization 
dynamics with the one of the transient charge asymme¬ 
try, which is extracted at selected time delays from the 
normalized intensity imbalance between the TSS bands 
at opposite +fc||,a and —fc||,a; wave vectors (see Supple¬ 
mentary). The asymmetry is evaluated at an energy of 
^0.3 eV, in correspondence with the spin-resolved mea¬ 
surements of Figs. 3a and 3b. We clearly observe that 
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FIG. 4: Schematics of the microscopic mechanisms underiying the uitrafast reiaxation of the Dirac-cone spin 
poiarization. From (a) to (d), subsequent stages of the spin-relaxation process are shown. The straight red lines and the 
parabolic black curve represent the TSS and BCB above the Fermi level, respectively. Blue filled circles denote photoexcited 
electrons, and green arrows correspond to the out-of-plane spin polarization, (a) A circular pump-pulse of positive helicity 
(C-b) creates a thermalized spin-polarized population of hot electrons in the TSS. (b) The initial relaxation of the TSS spin 
polarization mainly proceeds through intraband scattering processes, (c) Once electrons within the BCB decay to their band 
minimum, the unpolarized BCB electron population is transferred into the Dirac bands via phonon-assisted decay, so that 
bulk-to-surface coupling further contributes to the spin depolarization, (d, e) Subsequently, this process continues until the 
Dirac cone is largely populated by BCB electrons, leading to a suppression of the charge current on the surface (see text). 


the timescales associated to spin decoherence and relax¬ 
ation of the charge asymmetry within the Dirac cone are 
drastically different. In particular, up to ^600 fs, while 
the spin polarization is considerably reduced, the charge 
asymmetry remains nearly constant. This behavior re¬ 
veals that intraband scattering is a spin-depolarizing pro¬ 
cess that develops under the preservation of the tran¬ 
sient charge current on the surface. At the same time, 
the electron bottleneck effect arising at the Dirac node 
slows down the depolarization process while maintaining 
the charge asymmetry unchanged. After ^600 fs, once 
the BCB electrons start to more effectively decay into 
the TSS bands trough interband scattering, the spin po¬ 
larization is further reduced, and the charge asymmetry 
continuously decreases until it nearly disappears within 
~1.85 ps. Note that the decay time is slightly shorter 
than the ~2.5 ps needed for the TSS population to re¬ 
lax below ^0.3 eV (see Fig. If), suggesting that the 
quenching of the charge current occurs once the Dirac 
cone is largely filled by electrons originating from the 
BCB. From these results we derive that the slow charge- 
current decay is directly caused by the interplay between 
the electron bottleneck arising at the Dirac node and 
BCB-emptying processes, which simultaneously act as 
spin-depolarization channels. Hence, the suppression of 
the spin polarization is caused by the cooperative action 
of intraband and interband scattering processes. We em¬ 
phasize that in a nonequilibrium massive-fermion state, 
we would have expected the charge current to immedi¬ 
ately drop within the electron-electron scattering time 
due to elastic backscattering under 180 degrees, differ¬ 
ently from what is observed here. The slow charge dy¬ 
namics in Fig. 3c additionally suggests that a transient 
gap at the Dirac node does not open during consecu¬ 
tive momentum-relaxation stages. This fact is consistent 
with our general interpretation concerning the spectral 
weight suppression at the Dirac node. 

Based on this conclusion, in Fig. 4 we illustrate a 
simplified scheme that focuses only on the mechanisms 
underlying the spin-relaxation process. Up to ^600 fs, 


the oriented spin-polarized TSS population (Fig. 4a) de¬ 
cays mainly through intraband scattering (Fig. 4b). The 
electron momentum relaxation in the presence of spin- 
orbit coupling results in small random rotations of the 
electron spins due to scattering with phonons. In this 
process, the spin-orbit field can be viewed as an effec¬ 
tive magnetic field acting on the electron spins. Every 
scattering event results in a randomization of the elec¬ 
tron spin-precession axis around this field, causing spin 
depolarization. Subsequently, after the BCB electrons 
have reached their band minimum, the TSS spin polar¬ 
ization additionally decays through interband scattering 
(Fig. 4c). We emphasize that because the TSS and BCB 
edges are close to each other in energy-momentum space, 
interband scattering is likely to occur before the BCB 
electrons decay to their band minimum, but in a less ef¬ 
fective way. The process can be viewed as transfer of 
unpolarized electrons into the TSS bands, implying that 
bulk-to-surface coupling and intraband scattering coop¬ 
eratively depolarize the electron spectrum within ^1.2 
ps (Fig. 4d). The subsequent dynamics proceeds once 
the TSS population is completely filled with unpolarized 
BCB electrons at ~1.85 ps (Fig. 4e). These electrons 
remain unpolarized throughout the decay towards the 
Fermi level due to the strong contribution of unpolarized 
bulk-valence band states located right below and in the 
immediate vicinity of the Fermi level (see Fig. le). 

Our results establish that the transfer of BCB elec¬ 
trons into the TSS is a spin-depolarizing event, an effect 
that we attribute to the complex nature of the phonon- 
mediated interaction between the bulk and surface elec¬ 
tronic wave functions. The partial filling of the TSS 
bands with unpolarized electrons respects the total wave 
function symmetry of the TSS, as it is composed by 
different spin orientations on different atoms gOl SI]. 
Our photoemission calculations in equilibrium addition¬ 
ally reveal that the strength of the bulk-surface interac¬ 
tion is possibly mediated by the surface resonances ap¬ 
pearing at the BCB edges |43]. These observations re¬ 
call a strong need for refined spectroscopical calculations 



















of time-resolved photoemission spectra, which quantita¬ 
tively describe dynamical processes as a function of time. 
A first theoretical description of time-resolved photoe¬ 
mission for real systems was introduced recently |46| . 
Although the numerical implementation will need some 
time, at present such type of calculations are not pos¬ 
sible using only perturbative methods. It would be in¬ 
teresting to perform similar studies as the one presented 
here using lower pump photon energies so that all elec¬ 
trons are excited within the bulk band gap, or to use 
more bulk-insulating samples. Another alternative is to 
use ultrathin TI films without the contribution from the 
bulk conduction. 

To summarize, we have shown that circularly-polarized 
infrared pulses of opposite helicities allow all-optical 
control of the nonequilibrium spin polarization of TSSs 
on ultrafast timescales. The dynamics of the generated 
spin-polarized electrical currents exhibits two distinct 
timescales for the spin and charge degrees of freedom 
on the surface. We have provided key evidence for the 
mechanisms underlying this behavior. Our observations 
establish the utilization of photoexcited TSSs as unique 
channels in which to drive ultrafast currents with fully 
controllable spin configurations. The results presented 
here are of critical importance in understanding the spin 
properties of TIs following a perturbation by intense 
laser helds and open the way for novel applications of 
these materials in ultrafast spintronics. 

Methods 

Spin-resolved tr-ARPES experiments are performed at 
room temperature in ultrahigh vacuum better than 1 • 
10“^° mbar. We use a Mott-type spin polarimeter cou¬ 
pled to a Scienta R8000 hemispherical analyzer, allowing 
us to detect the out-of-plane and one tangential com¬ 
ponent of the spin polarization (see Supplementary for 
more details). As pump and probe pulses we use the 
hrst and fourth harmonics derived from a Coherent RegA 
Ti:Sapphire amplified laser operating at 150 kHz. Sb 2 Te 3 
single crystals are grown by the Bridgman method and 
cleaved in situ. Photoemission calculations are based on 
multiple scattering theory within the one-step model of 
photoemission including wave-vector, spin and energy- 
dependent transition matrix elements mi SSI. We use 
a fully-relativistic version that is part of the SPR-KKR 
program package [49l [50], with spin-orbit coupling in¬ 
cluded self-consistently. 
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